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Abstract. Strong modification of event-by-event fluctuations in global observables was proposed by several 
authors as an indicator of the phase transition in system evolution. We report the results of the analysis of 
STAR TPC data on the charged multiplicity and mean transverse momentum fluctuations in Au+Au collisions at 
\Aaw = 130 GeV. We discuss the centrality dependence of the observed fluctuations and compare the results with 
available data from lower collision energies and ISR data. 



The question of the hadronization of the system created in ultrarelativistic heavy ion collisions is one of the most 
interesting questions in this field. Does the system evolution include the phase transition? Event-by-event fluctuations 
of the mean transverse momentum and/or multiplicity, are considered to be one of the important tools to identify such a 
phase transition. The fluctuations depend on the nature of the phase transition. A first order phase transition may lead to 
large fluctuations in energy density due to formation of QGP droplets[l, 2]. Second order phase transitions may lead to 
divergence in the specific heat; it would also increase the fluctuations in energy density due to long range correlations 
in the system [3, 4]. One could observe them as fluctuations in mean transverse momentum if matter freezes out at 
the critical temperature T c [5, 6]. The system freeze-out near the QCD tri-critical and critical points would also lead to 
change in the fluctuations pattern [3, 4, 7]. There exist also other sources of the abnormal dynamical fluctuations. The 
formation of the disoriented chiral condensate [8] should lead to large fluctuations in the ratio of neutral to charged 
pion yields at low p,, and consequently to fluctuations in (p t ) of charged particles. Jets may give rise to event-by-event 
fluctuations in the high p t region. 

For the multiplicity fluctuations, predictions range from enhanced multiplicity fluctuations connected to the produc- 
tion of QGP droplets and nucleation processes in a first order phase transition, to a strong suppression of fluctuations as 
a consequence of rapid freeze-out just after the phase transition. In this case, final state values of conserved quantities, 
such as net electric charge, baryon number, and strangeness would not be strongly modified from their value in the 
QGP stage. The fluctuation of these quantities could be reduced by a factor 2 - 4 if a QGP is produced [9, 10]. The 
production and size of QGP droplets may critically depend on the collision impact parameter, leading to centrality 
dependence of the fluctuations to be an important observable. If the fluctuations are due to the particle production via 
any kind of clusters (e.g., resonances, strings (mini)jets, independent AW-collisions, etc.) and the relative production 
of such clusters do not change with centrality, the correlation strength should be inversely proportional to the particle 
multiplicity. New physics should appear as a deviation from such a dependence. 

Here, we report on our measurement of charged particle multiplicity and mean p t fluctuations as a function of 
collision centrality in Au + Au collisions at an energy of ^/s^n = 130 GeV. Specifically, for the case of multiplicity 
fluctuations we discuss fluctuations in the difference of the relative number of produced positive and negative charged 
particles measured in a fixed rapidity range: 



INTRODUCTION 




(1) 



where N a and Nb are the corresponding multiplicities. One can consider fluctuations in the production of charged 
particles, N+ and N-, as well as in the specific cases of proton and anti -proton, N p and N-p, and positive and negative 



kaons, Nk+ and Nk- ■ The former amounts to a measurement of net electrical charge fluctuations, whereas the latter 
cases are related to net baryon number and net strangeness fluctuations, respectively. The quantity v a b used in our 
analysis can be related to the fluctuations in multiplicity ratios [12, 11], but have the advantage of not being limited to 
large particle multiplicities. 



Statistical and Dynamical fluctuations 

The fluctuations of any observable measured on an event-by-event basis consists of two parts: statistical fluctuations, 
the uncertainty in the measurement of a given quantity due to finite event multiplicity, and dynamical fluctuations due 
to correlated particle production, the subject of interest. The statistical fluctuations are the ones one would observe in 
the case of independent particle production provided that the single particle inclusive quantities remain the same as in 
real events. 

Multiplicity fluctuations. A difficulty inherent in the interpretation of multiplicity fluctuations is the elimination 
of effects associated with uncertainties in the collision centrality, often referred to as volume fluctuations. They, 
in particular, induce a positive correlation in particle production which does not depend on the intrinsic dynamical 
properties of the colliding system. Fluctuations in the difference of relative multiplicities, v a b, defined in Eq. 1, are 
however free from this problem. This analysis is thus restricted to the study of such relative multiplicities. 

The magnitude of the variance, V a b, is determined by both statistical and dynamical fluctuations. Taking into account 
the fact that for independent particle production the multiplicity fluctuations follow the Poisson distribution, the 
dynamical fluctuations can be evaluated as 

Vab.dyn = V 'ab — Vab.stati ^ab.stat = , + , • (2) 

(Na) (N h ) 

The value of V a b,dyn is determined by intrinsic correlations between produced particles averaged over the pseudo- 
rapidity and transverse momentum range under study. It is, to first order, independent of the detector acceptance 
(provided it is smaller than typical long-range correlation scale of 1 - 2 units of rapidity) and particle detection 
efficiency as can be seen from its relation to the inclusive two and single particle densities: 

Vab,dyn = Raa + Rbb — 2R a b, (3) 

where R a b is a correlation function similar to that used in multi-particle production analysis [13, 14, 15] and defined 

as 

R b = .LYP2,ab{T\a,T\b)dT\gdT\ h _ ^ 
lAYPlA^a)dy\a.f A YPl,b(^b)dy\b 

where Pi(r|) = dn/dr\, and P2(ila, r lfo) = d 2 n / dr\ a dr\b are single and two particle (pseudo)rapidity densities, respec- 
tively. 

Mean transverse momentum fluctuations. We use the following notation for the mean transverse momentum in 
each event and average over all events in a given event sample (centrality bin): 

1 M i N ev , 

(P')=mLp^ <fo» = irl>*>y (5) 

The fluctuations in {p t ) becomes 

= ((^> 2 >-((A>> 2 -i((Eft 2 1 + E^r, / ))-((A)) 2 (6) 

a 2 + M ~ 1 a 2 with a 2 - {{P ' }) - {{pt))2 - (T) 

- G{p,),stat + —jfi- a {p,),dynam' wltn °(p,).stat ~ - ^ ■ V) 

where M is the multiplicity. In Eq. 7 we have introduced the notation for the non-statistical (dynamical) fluctuations 

in (p t ) as: 

a }p,),dynam = ((Pt,iPt,j)ift) - {{Pt)Y \ (8) 

which is zero in the case of no correlation between particle transverse momenta. As above, in this equation the double 
angle brackets represent the average over particle pairs in a specified (pseudo)rapidity and transverse momentum 
window, and the average over all events with a given multiplicity. 
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FIGURE 1. (left) v + _ </„, as a function of centrality (event raw multiplicity), and (right) v+_/v + _ 5 , fl , ~ 1 + {{N+) + 
(yv_))v + _. rfv „/4. 



RESULTS 

STAR detector and data selection. The STAR detector consists of many detector subsystems in a large solenoidal 
magnet with the Time Projection Chamber (TPC) as the main tracking device. Details about the STAR TPC can be 
found in reference [16]. In the year 2000, the TPC was operated in a 0.25 T magnetic field allowing tracking of particles 
with p t greater than 75 MeV/c. The experimental setup includes also the scintillator trigger barrel (CTB) surrounding 
the TPC and two Zero Degree Calorimeters (ZDCs) located at 18 m from the center of the TPC. The ZDCs measure 
mainly fragmentation neutrons and are used in coincidence as the experimental trigger for minimum bias events. 

We determine the centrality of the collision by the total reconstructed charged track multiplicity, N c h, in the 
pseudorapidity region |r|| < 0.75. Eight centrality bins used in this analysis are the same as used in [17] and correspond 
to the fraction of the total geometrical cross section of the top 6%, 11%, 18%, 26%, 34%, 45%, 58%, and 85%, 
respectively. 

A total of about 140 thousand minimum bias events with a primary vertex position within 75 cm of the center of the 
TPC have been used in this analysis. The tracks are required to pass within 3 cm of the primary vertex and the ratio of 
the number of points on a track to the maximum possible number of points is required to be greater than 0.5 in order 
to avoid double counting of tracks due to track splitting. In this analysis we use only tracks within the pseudo-rapidity 
range —0.75 < r| < 0.75 and transverse momentum 0. 1 < p t < 2.0 GeV/c, where the track reconstruction efficiency is 
estimated to be 90 ±10%. 

Multiplicity fluctuations. Fig. 1 shows the dynamical fluctuations, V+_, ( /->,„, in the pseudorapidity region —0.5 < r| < 
0.5, as a function of the total multiplicity N c h- The dynamical fluctuations are finite and negative: a clear indication 
of the correlation in the production of the positive and negative particles. We verified that, within statistical errors, 
the same value is obtained when the pseudorapidity regions used to count positive an negative tracks were separated 
by a Ar| = 0.125 gap, which should eliminate problems related to track splitting and resolution of two close tracks. 
We observe less than 10% decrease in the strength of the correlation when varying the pseudo-rapidity range from 
|T|| < 0.125 to |T|| < 0.75. 

The absolute strength of the dynamical fluctuations decreases monotonically with increasing collision centrality. 
In order to check the deviation from the inverse proportionality to multiplicity, in Fig. lb we present the ratio of the 
observed fluctuations to the statistical fluctuations'. This ratio is related to the quantities cog and D used in other 
works [9, 10, 18] as \ a bftab,stat ~ <*>£> = &Q 2 /N c h = 0.25D. The relationship is strictly exact if (N+) = (N-). The ratio 
has a value of about 0.8 nearly independent of the collision centrality in qualitative agreement with the expected value 
0.75 based on a thermalized resonance gas as calculated in [11] for particles produced in a rapidity interval Ay = 1. 
The measurement, however, deviates drastically from calculations based on a rapid QGP freeze out [10, 9], which 
predict values of «0. 2-0.4. 



Note, that this ratio depends on the size of the acceptance used for the measurement. 



We consider these results in the light of correlation functions measured in pp collisions at CERN [13, 20, 21] 
with the use of Eq.3. We account for the unavailability of measurements at the same energy by interpolating the data 
obtained at lower and higher collision energies (parameterization from [22]). We also note that the correlation function 
for opposite charge particles, R+-, is found to be approximately twice as strong as the same sign particles correlations, 

R ++ pa R [13, 14], and independent of the collision energy. Using the single charge particle pseudorapidity density 

pi(r| = 0) pa 2.05 and C2(0,0) = p2(T|i = 0,r|2 = 0) - pi(r|i = 0)p i (t|2 = 0) pa 2.8 one finds charged-charged 
correlation function R cc pa 0.66. Furthermore, assuming equal multiplicities of positively and negatively charged 
particles, one finds R cc pa (R++ +/?__ +2/? + _)/4 pa 1.5R++, which we use to estimate the correlation measured 
in this work as R ++ +R-- — 2R + - pa — 2R ++ pa 4R cc /3 pa 0.88. Under the assumption that the correlations are due 
to production in a finite number of sources (clusters), the correlation function should be inversely proportional to the 
particle density. In central Au+Au collisions, the charged particle density is about 550 compared to approximately 2.06 
in pp collision. Such a dilution would give for the correlation function a value of 0.88 • 2.06/550 pa 0.003, a number 
which is about 30% larger than we observe (see Fig. la). 

Mean p t fluctuations. The results for the mean transverse momentum fluctuation in the region —0.75 < r| < 0.75 
are presented in Fig. 2. Selecting specific pairs for the averaging in Eq. 8 one can suppress or enhance the contribution 
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FIGURE 2. 

multiplicity. 

to (p t ) fluctuations of different origin. We have also correlated: (1) particles in pseudo-rapidity region —0.75 < r| < 
—0.05 with particles in 0.05 < r| < 0.75. The 'gap' between the two regions eliminates effects such as Bose-Einstein 
correlations (HBT) or Coulomb final state interactions. (2) Positive particles with negative particles, which is expected 
to enhance the contribution from resonances. No significant difference has been observed. The somewhat smaller value 
(less than 10%) in the case of the correlation of particles separated by the gap in pseudorapidity could be explained 
by larger average distance in the pseudorapidity space between the two particles in a pair. The relative fluctuation 
in mean transverse momentum was found to be a dynam I ((Pt)) ~ 1.2 ±0.2%. The (relative) systematic error in our 
measurements of o? > dynam * s arx>ut 5%- It nas been estimated by varying track selection cuts, primary vertex position 
within the TPC, and comparing the Monte-Carlo results before and after processing through the detector simulation. 

Event-by-event dynamical fluctuations have also been analyzed by several experiments using the so called <P Pt [23] 
measure (the approximate relation to of i dynam * s ta k en from [19]): 



^({( Pt )-Mplf)/(M)- 



* pt, inclusive ' 



(pt) ,dynan. 



2a 



Pi, inclusive 



(9) 



and close to it the difference factor Ao P; [24]: 



' pt , inclusive ' 



(10) 



For the 6% most central collisions STAR measures [24] Aa Pt pa 35 MeV, which is consistent with results reported 
here for (S^^ynam ( one can check this with the help of an approximate expression given by Eq. 9 and using the 
other parameters presented in [24]). In Pb+Pb and Pb+Au collisions at CERN SPS (^/snn — 17 GeV) <P Pt have been 



measured by NA49 [25] and NA45/CERES [26]. NA49 reported [25] <$> Pt = 0.6 ± 1 MeV/c for the rapidity region 
4.0 < y < 5.5 and 5% most central collisions. The observed fluctuations are extremely small, but it would be incorrect 
to compare these numbers with the STAR measurement, since NA49 results were obtained in the forward rapidity 
region, while STAR measurements are at midrapidity. The CERES collaboration has measured the fluctuations in the 
central rapidity region. They report <P Pt = 7.8 ±0.9 MeV/c in central Pb+Au collisions [26]. As seen from Eq. 9, 
Op, is directly proportional to the number of reconstructed tracks used for its calculation (subject to acceptance 
cuts and tracking efficiency) which complicates the comparison of the results from different experiments. For a 
rough comparison one can take into account the fact that the CERES multiplicity ~130, ((pt)) «420 GeV/c and 
<V, inclusive ~ 0.270 GeV/c. Then <P Pl w 8 MeV corresponds to o^ pt)4ynam / {{p,)) w 1 .4%. 

The dynamical part of mean transverse momentum fluctuations has been measured at the ISR [27] in pp collisions. 
This was done by analyzing the multiplicity dependence of G{ Pt ) under the assumption that in pp collisions the 
dynamical part in (p t ) fluctuations does not depend on multiplicity. It was observed that <5(p t ) 4yna ml {(Pt)} ~ 12%. 
Rescaling of this quantity with (the square root of) the ratio of multiplicity densities in pp and Au+Au collisions 
yields the fluctuations in Au+Au: <5^ 4ynam l {{pt)) ~ 0.8%, about 50% less than observed in this experiment. 

To understand the role of (mini)jet production in (p t ) fluctuations we have performed an analysis of central HIJING 
events with hard processes switched "on" or "off" and with the same cuts as in the real experimental data. The 
jets increase the fluctuations from 0( Pt },d yna m/ ((Pt)) ~ 1.0% to about 1.1%. Note that both numbers are very close 
to the experimentally observed value. Contrary to HIJING, the RQMD event generator includes secondary particle 
rescatterings, but yields approximately the same values as HIJING (within 20%, limited by our statistics of the RQMD 
events). 

In summary, we have reported measurements of the multiplicity and mean transverse momentum fluctuations 
in Au+Au collisions at RHIC. While a clear correlation in particle production has been measured, no qualitative 
difference from the lower energy AA collision experiments as well as pp collisions at similar energies has been 
observed. The centrality dependence of the observed correlations is consistent with being inversely proportional to 
total event multiplicity, the dependence expected in particle production via clusters. 
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